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Abstract 

We study the semi-discrete directed random polymer model introduced by O’Connell 
and Yor. We obtain a representation for the moment generating function of the 
polymer partition function in terms of a determinantal measure. This measure is 
an extension of the probability measure of the eigenvalues for the Gaussian Unitary 
Ensemble (GUE) in random matrix theory. To establish the relation, we introduce 
another determinantal measure on larger degrees of freedom and consider its few 
properties, from which the representation above follows immediately. 


1 Introduction 

In this paper we consider a directed random polymer model in random media in two (one 
discrete and one continuous) dimension introduced by O’Connell and Yor [59]. For N 
independent one-dimensional standard Brownian motions Bj(t ), j = 1, • • • ,N and the 
parameter /?(> 0) representing the inverse temperature, the polymer partition function is 
defined by 

Z N (t)= f e P( B l(s 1 )+ B 2{s 1 ,S 2 )+-+ B N(s N - 1 ,t)) dsi . . . dsjv-l- (1.1) 

J 0<si <-<Sff-i<t 

Here Bj(s,t ) = Bj(t ) — Bj(s ), j — 2, • • • , N for s < t and — Bi(si) — B 2 (si,s 2 ) — • • • — 
Bn(sn-i, t) represents the energy of the polymer. In the last fifteen years much progress 
has been made on this O’Connell-Yor polymer model, by which we can access some explicit 
information about Z^{t) and the polymer free energy F/v(t) = — log(Zjv(f))//l [71110111111361 
SHS5HI71I52E71E2]. The first breakthrough was made in the zero temperature (/3 -A oo) 
case. In this limit, — F/v(f) becomes 

/jv( 1) := - lim F N (t) = max (Hi(si) + B 2 (s 1 ,s 2 ) 4-b B N (s N -i,t)) (1.2) 

p—> OO 0<Sl<-”<SjV-l<£ 
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where — /at(£) is the ground state energy. For /w(t), the following relation was estab¬ 
lished 0136]: 


Prob (/at(X) < s) 


' (—oo,s] N 


N 

n dx i 

3 = 1 


• Pgue(^i, • • • , t ), 


N 

Pgue(^i, ■ ■ • ,X N -,t) = P] 

3 = 1 


e~ x V 2t 

j\P~'V2rt 


n ( Xk - 


(1.3) 

(1.4) 


where -Pgue(Ti, • • • , x N]t) is the probability density function of the eigenvalues in the 
Gaussian Unitary Ensemble (GUE) in random matrix theory [31 1351152] . This type of 
connection of the ground state energy of a directed polymer in random media with random 
matrix theory was first obtained for a directed random polymer model on a discrete space 
Z+ [H] by using the Robinson-Schensted-Knuth(RSK) correspondence. Eq. (II.3ft can be 
regarded as its continuous analogue. Note that (II.41) is written in the form of a product of 
the Vandermonde determinant rii<j<A:<Ar( x fe — x j)- This feature implies that the m-point 
correlation function is described by an m x m determinant, i.e. the eigenvalues of the GUE 
are a typical example of the determinantal point processes 0- In addition based on this 
fact and explicit expression of the correlation kernel, we can study the asymptotic behavior 
of / 7 v(f) in the limit N —> oo. In 036], it has been shown that under a proper scaling, the 
limiting distribution of /jv(i) becomes the GUE Tracy-Widom distribution [ 75] . 

In this paper, we provide a representation for a moment generating function of the 
polymer partition function (11.11) which holds for arbitrary P(> 0): 


E 


exp 


•- pu z N (t) 

p2{N—l) 


N 


Y[ dx j fF{ x i ~ U) ■ W(xi, ■ ■ ■ ,X N ]t), 


3=1 


N 


W{x-i, ■ ■ ■ ,x N -t) = (x k - Xj) ■ det (il>k-i(xj',t ))^ k=1 , 

j=l^' 1 <j<k<N 

where /f(^) — l/(e^ + 1) is the Fermi distribution function and 


1 P°° 

M x ;t) = ^ dwe - iwx - wH / 2 


\IW) 


r (i + iw/p) 


N ■ 


(1.5) 

( 1 . 6 ) 


(1.7) 


For more details see Definition |Tj and Theorem [2] below. This is a simple generalization 
of (11.31) to the case of finite temperature. We easily find that it recovers (11.31) in the zero- 
temperature limit (P —>■ oo). Note that the function W{x\, ■ ■ ■ ,x ^; t) is also written as a 
product of two determinants and thus retains the determinantal structure in (ll.4|) . 

In most cases, to find a finite temperature generalization of results for zero-temperature 
case is highly nontrivial and in fact often impossible. But for the O’Connell-Yor polymer 
model and a few related models, rich mathematical structures have been discovered for 
finite temperature and the studies on this topic entered a new stage [2l ll0ll251[37ll57ll61I[66 - 
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[69]. In [57], O’Connell found a connection to the quantum Toda lattice, and based on the 
developments in its studies and the geometric RSK correspondence, it was revealed that 
the law of the free energy F^(t) is expressed as 


Prob (—F N (t ) + N 1 log/3 2 < s\ = [ clx± [ TT dxj ■ m(x 1 , ■ ■ - , x N -, t). (1.8) 
VP/ J(-oo,s] Jrn-i 

Here the probability measure m(x 1 , • • • , xn] t ) riy=i ^Xj, which is called the Whittaker mea¬ 
sure, is defined by the density function m(x 1 , • • • , xn ; t) in terms of the Whittaker function 
• • • ,xn) (for the definition, see [5T]J) and the Sklyanin measure SN(X)dX (see (12. lOh 
below) as follows, 

m(x 1 , • • • , x N ] t) = ^ 0 (^X 1 , • • • , (3x N ) / dX if!_x/p(Px 1 , • • • , px N )e^ =1 ^ t/2 s N ( A//3), 

J(i R) n 

(1.9) 

where A represents (Ai, - - - ,Xn)- In contrast to (11.41) . the density function (11 .9(1 is not 
known to be expressed as a product of determinants and the process associated with (11.9(1 
does not seem to be determinantal. Nevertheless some determinantal formulas for the 
O’Connell-Yor polymer have been found: First in [57], O’Connell showed a determinantal 
representation for the moment generating function (LHS of (11.5ft ) in terms of the Sklyanin 
measure. (See (I2.9[) below.) Next in [10], Borodin and Corwin obtained a Fredholm de¬ 
terminant representation for the same moment generating function (see (14.231) below). A 
direct proof of the equivalence between the two determinantal expressions was given in [13] . 
In pro], by considering its continuous limit, the authors also obtained an explicit represen¬ 
tation of the free energy distribution for the directed random polymer in two continuous 
dimension described by stochastic heat equation (SHE) [TUI ITT] , The distribution in this 
limit, which describes the universal crossover between the Kardar-Parisi-Zhang (KPZ) and 
the Edwards-Wilkinson universality class, was first obtained in [211661 169] and can be inter¬ 
preted also as the height distribution for the KPZ equation [44] , Furthermore in ra. they 
consider not only the O’Connell-Yor model but a class of stochastic processes having the 
similar Fredholm determinant expressions, the Macdonald processes, the probability mea¬ 
sures on a sequence of partitions which are written in terms of the Macdonald symmetric 
functions and include the Whittaker measure defined by (11.81) as a limiting case. 

The purpose of this paper is to investigate further the mechanism of appearance of 
such determinantal structures and (ll.5|) is the central formula in our study. Although 
W(x 1 , • • • ,x N ;t) nJLi dxj defined by (11.61) is not a probability measure but a signed mea¬ 
sure except when /3 —» 00 , a remarkable feature of this measure is that it is determinantal 
for arbitrary f3 contrary to the Whittaker measure (11.91) . This determinantal structure 
allows us to use the conventional techniques developed in the random matrix theory and 
thus from the relation we readily get a Fredholm determinant representation with a kernel 
using biorthogonal functions which is regarded as a generalization of the kernel with the 
Hermite polynomials for the GUE. In (11.61) . the parameter /3, which originally represents 
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the inverse temperature in the polymer model appears in the Fermi distribution function 
fp{x — u) with the chemical potential u as well as i])k(x\t) (11.7ft in RHS. This fact with 
the determinantal structure suggests that the RHS might have something to do with the 
free Fermions at a finite temperature. Related to this, a curious relation of the height of 
the KPZ equation with Fermions has been discussed in [28] . 

For establishing the relation, we introduce a measure on a larger space R JV < A, + 1 )/ 2 . By 
integrating the measure in two different ways, we get its two marginal weights. In one 
formula appears a determinant which solves the N dimensional diffusion equation with 
some condition (see (12.111) . (13.61) . and (13.71) ) and the other one with a symmetrization is 
exactly the RHS of (1 1.5 [) . The relation (1 1.5 j) follows immediately from the equivalence of 
these two expressions. Our approach is similar to the one by Warren [78] for getting the 
relation (II.3ft . Actually in the zero-temperature limit /3 —> oo, we see that the integration 
of the measure is written in terms of the probability measure introduced in [78], which 
describes the positions of the reflected Brownian particles on the Gelfand-Tsetlin cone. 
Note that the Macdonald processes (especially the Whittaker process in our case) [TO] are 
also another generalizations of [78] . Although the Whittaker process has rich integrable 
properties, they do not inherit the determinantal structure of [78]. On the other hand, our 
measure is described without using the Whittaker functions and keeps the determinantal 
structure. Furthermore combining (11.51) with the fact that the quantity can be rewritten 
as the Fredholm determinant found in [10] (Corollary fl3l and Proposition fl5l below), our 
approach can be considered as another proof of the equivalence between (14.231) and (12. 9j) 
hi [13]. One feature of our proof is to bring to light the larger determinantal structure 
behind the two relations. 

This paper is organized as follows. In the next section, after stating the definition of 
a determinantal measure, we give our main result, Theorem [7] and its proof. The proof 
consists of two major steps: we first introduce in Lemmata determinantal representation 
for the moment generating function which is a deformed version of the representation (12.9[) 
in [57]. Next we introduce another determinantal measure on larger space R iV ( JV+1 )/ 2 and 
then we find two relations about its integrations which play a key role in deriving our main 
result. In Sec. 3 we show that this approach can be considered as an extension of the 
one in [78]. In Sec.4., we consider the Fredholm determinant formula with biorthogonal 
kernel obtained by applying conventional random matrix techniques to our main result. 
The scaling limit to the KPZ equation is discussed in Sec.5. We check that our kernel goes 
to the one obtained in the studies of the KPZ equation. A concluding remark is given in 
the last section. 


2 Main result 

In this section, we introduce a measure W(x i, • ■ • , x^', t) IdyLi ^ x j CPU) , state our main 
result and give its proof. 
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2.1 Definition and result 

Definition 1 . Let ifk(x', t), k — 1,2, be 


k(x]t) = 


2tt 


dwe~ iwx ~ w2t/2 - 


(iw 


r (i + iw/p) 


N ' 


( 2 . 1 ) 


For (xi, • • • ,x N ) e 


pN 


a function W(xi, ■■ ■ ,xjsr;t ) is defined by 

N „ 


,x N ;t) ’ ]^[ (x m - xf) ■ det {i>j-i{x k \t))^ k=1 . 


j =1 J 1 <l<m<N 


( 2 . 2 ) 


Remark. We find that - ,xw]t) is a real function on ~R N , since by definition 

t) is real for any k = 0,1, 2, • • • , N — 1, P > 0 and t > 0. But in general, the positivity 
of this measure is not guaranteed. For example ifo(x',t ) shows a damped oscillation and 
can take a negative value for some x. Thus at least for the case N — 1, W(x,t ) = i])q{x\ t ) 
can be negative. 

We discuss the zero-temperature limit /3 —* oo of hF(xi, • • • ,x^\t). Noting T(l) = 1, 
we see 

i r°° p—x 2 /2t / -| \ § 

= a? = 7ST (s) ft 



where we used the integral representations of the nth order Hermite polynomial H n (x) (see 
e.g. Section 6.1 in 0). 


H n (x) = ( 2 f" f 
J— oo 


(2.4) 


Note that {t/2) k ^ 2 Hk {x/\pl it) is a monic polynomial (i.e. the coefficient of the highest 
degree is 1) and 


^ e ~ x }l 2t 

lim det (^-lfo; t))f fc _ x = TT . - • TT (x fc - xf). (2.5) 

0 -+°° 3 ’ . "7 t- 7_i v / 2vrt AA 

J=1 v 1 <j<k<N 

Thus we find 

lim W(xi, ■ ■ ■ ,x N ;t) = P G ue(^i, • ■ ■ ,x N ]t), (2.6) 

/3—too 

where Tgue^i, • • ■ ,X]y',t) is defined by (II.4jl . The function W{x i, ■ ■ • ,xjsr]t) can be re¬ 
garded as a deformation of (11.41) which keeps its determinantal structure. 

In this paper, we provide a determinantal representation for the moment generating 
function of the polymer partition function (11.11) in terms of the function (12.21) . 
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Theorem 2. 


E [e • )2l ' v J = J N \\dxj f F (xj - u) ■ W(xi, ■ ■ ■ , x N ; t) (2.7) 

where f F (x) = l/(e^ x + 1) is the Fermi distribution function. 

By (II.2p . (12.6p and the simple facts 

lim e~ ePx = lim f F (x ) = ©(— x), (2.8) 

/3—>oo /3—yx 

we find that the zero temperature limit of (12.71) becomes (11.31) . 

Because of the determinantal structure of W{x\, ■ ■ ■ ,x^;t), we can get the Fredholm 
determinant representation for the moment generating function by using the techniques in 
random matrix theory. Recently another Fredholm determinant representation has been 
given based on properties of Macdonald difference operators [10]. The equivalence between 
them will be shown in Sec. [H 


2.2 Proof 


Here we provide a proof of Theorem [2l Our starting point is the representation for the 
moment generating function given in [ 57]: 


E 



Pu Z N (t) 
pHN-1) 




u\j+\?t/2 p 



where 0 < e < (3 and Sjv(A)dA is the Sklyanin measure defined by 


sn(X) 


1 

(2 7ri) N N\ 


n 


sin 7 t(A i 

71 


— n _ • 

i>j 


(2.9) 


( 2 . 10 ) 


This relation was obtained by using the properties of the Whittaker functions [22l[74] and 
the Whittaker measure (14.211) . 


Lemma 3. 


e~^Z N (t) 

E e p 2(N ~ 1] = 


N 

/ T\dx e f F (xt-u) -G(x ir -- ,x N ]t) 
Jr n £=1 


where f F (x) is defined below (14 .1 7j) and 

G{x i, • • • , x N ] t) = det (F jk (x N - j+ p t))f >k=1 , 

Fjk(x;t ) = 


dX e ~ Xx+xH / 2 /tt 7tAV _1 

-»U cot y) 


2vr ir, i x 


F(f + 1 


with 0 < e < /3. 


( 2 . 11 ) 


( 2 . 12 ) 

(2.13) 
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We will discuss an interpretation of (12.12)1 in the next section. In this definition, we 
have arranged x^’s in the reversed order so as to relate (13.1711 . the zero-temperature limit 
of (12.121) . to the stochastic processes defined later in (13.2011 . 

Proof. Noting the relation 


j[ sin (ay — Xj) — sinay sin x 3 (cot x 3 — cot ay) 

l<i<j<N 1 


N 


N 


sin A 1 Xj ■ (cot ay — cot Xk) = siir v 1 Xj ■ det (cot f 1 ay) fc ( _ 1 , (2.14) 

3 = 1 


j=l l<k<t<N 

we rewrite RHS of (12.9)1 as 




N 


' SNl j 


1 

m 


N 


dX 


l-e) N J = 1 


n __j_ p—u\j + t/2-p / 

2m \ R 


A., 


P 


N 


7r \ \ N—l 


sin 


7r 


x det | ( — cot 77 A,- 

,/3 P 


k-\ 


N 


^ (A/' 1 )"., 


= det 


dX 

l-e 27TZ/3 f 


3 —«A+A 2 i/2-p / _ _ 

' /3 


j,fc=i 
N 


X\ N /onfA 


7T \ \ At—1 


7T 


i-i 


i cot i Al ^ 


At 




(2.15) 


where in the last equality, we used the Andreief identity (also known as the Cauchy-Binet 
identity) [4]: For the functions gj(x), hj(x), j = 1,2, such that all integrations 

below are well-defined, we have 


TFT/ f[ d Xj-det (g k (xj))" det {h k (xj))*= det ( [ dxgj(x)h k (x)) . (2.16) 

\Jr J j,k=l 


We notice that the factor e “ A r(— X//3) N (sm(TrX//3 )/tt) n 'in (I2.15P can be written as 


e -“ A r [ - 


A\ N ( s in f A 

P 


7r \ \ At—1 

A 


7T 


-1) 


At i Tj-g uX 


r i + 


aA — sin j|A 


(- 1 ) 


AT—1 r o° 


n —x\ 


r i +a 


At 


/3 


e P{x-u) _j_ 


dx 


(2.17) 


where we used the reflection formula for the Gamma function and the relation (14.31 jl . 
From (12.15)1 and (12.1 711 . we arrive at the desired expression (12.lip . □ 
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From (12. lip , we see that for the derivation of our main result (12.71) . it is sufficient to 
prove the relation 


r n r N 

/ T\dx t f F {xt - u) ■ G(xi, ■ ■ ■ ,x N ;t) = / TT dxj f F (xj - u) ■ W(x ll • • ■ , x N ; t). 
Jr n Jr» 


(2.18) 


where f F (x) is defined below (12.71) and W(x i, • • • , xjy] t ) is given in Deffiiition[H Note that 
this is a relation for the integrated values on M A '\ To establish this we introduce a measure 
on the larger space M. N ^ N+1 ^ 2 . 


Definition 4. Let x k be an array (x^\--- ,x^) where x^ - 
dx k = :1 dx^K We define a measure R u (x N ;t)dx N by 

R u (x n ; t) = JJ fi(x[ j) - 4^ 1} ) • det (f u {x 

1 <i<j<N 




G and 


(2.19) 


Here Xq 1j = u, Fij(x;t) is given by F t j(x: t) (j2T3j) with i = 1 and fi (x),i = 1,2, ••• 
is defined by using the Fermi and Bose distribution functions, /f(x) := \/[e^ x + 1) and 
fs{x) := l/(e^ x — 1) respectively as follows. 


fi(z) 


f F (x), i = 1, 
Ib(x), i > 2. 


( 2 . 20 ) 


Remark. The reason why both the Bose and Fermi distributions appear in our approach is 
not clear. The interrelations between them (see (12.281) - (12.30^ below) will play an important 
role in the following discussions. 

As in Fig [0 we usually represent the array x N graphically in the triangular shape. 
Although no ordering is imposed on x N , in the zero-temperature limit, R u (x N ]t ) has the 
support on the ordered arrays as in Fig. |T] (a) (see (j3.34p ). Fig. |Tj (b) represents the other 
ordered array called the Gelfand-Tsetlin pattern (see (I3.23P ). 


r 0) 

■' i 


r (3) 


X 


(3) 


7/ 


-7/ 


x 


( 2 ) 


X. 


( 2 ) 


7/ 


x 


(i) 


(b) 


r (3) 
1 1 


t (3) 

■ l 2 


X 


(3) 


R/ R/ 

t ( 2 ) t ( 2 ) 

A H/ 


x 


(i) 

i 


Figure 1: Triangular arrays (k = 3) (a) an element of 14 (13.341) (b) The Gelfand-Tsetrlin 
pattern (an element of (13.231) 1 




















As discussed later we will find that the moment generating function of the O’Connell- 
Yor polymer model is expressed as the integration of this measure R u (x N ; t ) over R N ( N+1 )/ 2 . 
We have other choices for the definition of R u (x N ; t) which give the same integration value. 
One example is 


Ru(x N ] 0 = 4 det - X i- 1 1] 

t=\ l ' 


l 

■ det 

i,j =1 



( 2 . 21 ) 


This comes form the following consideration. Let f sy m(x N ) be a function which is sym¬ 
metric under permutations of xf \ • • • , x'p for each j G {1, 2, • • • , N}. Then we see that 
Ru(x N ]t) (12.19P and R u (x N ',t ) have the same integration value: 


l«.N(N+l)/2 


dx_]y fsym(x_ n') Ru(x. JV) 


l R N(N+l)/2 


dx_]y fsym^^Liy') Ru(.X_ JV) 0* 


( 2 . 22 ) 


It can be shown as follows. From the symmetry of f S ym(x N ), LHS of the equation above 
becomes 


J R N(N+l)/2 

Here R u (x N : t) is defined by 


dx_]\f /sym(2ijv) R'u(.X.Ni 


J W N(N+l)/2 


dx_Nf 'sym(^ijv) Ru(x_N > ^)" 


(2.23) 


JV 


Ru{x n ; 2) = ^ 

£=i ‘ a u)eSj, j=i,-,N 


Ru {x N , t) , 


(2.24) 


where Sj is the permutation of 1,2, • • • , j and gf N denotes (x aW , • • • ,x a{N) ) with x a(3) = 
( X aU> n)>''' ’ We easily find the equivalence R u (x N \t ) = R u (x N ', t). Note that 


~0) _ -O'-i) 


> (i _ 1) ).det(F 1 .;(hS,y ) ;i))' V j _ i 


*) = n /i 

1 <i<j<N 
N j 

= SgncT (JV) n n fi (^0)(0 - • det i)). =1 


JV 

n 

i=i 


j=i *=i 
j 

- 0 ) 


= 11 sgnr'-- JJ./: 

2=1 


- 0 ) 

* rrW(i) 


jc-li 1 ’) • det (Fij(x) iV ' l ;t) 


„(Y.. 


JV 


m=i 


(2.25) 


Here in the last equality, G Sj, j = 1,2, , AT is defined by using cr (j b and cr01 

as crO _1 )r^)(A:) = a^(k), k — 1, • ■ • , j, where we regard cr0~b as an element of Sj with 
crO _1 )(j) = j. Further in the last equality we used u (iV ) = fl ^ =1 rO). Substituting (12.25]) 
into (12.2411 and using the definition of the determinant, we have R u (x N ;t) = R u (x N ;t). 

The function R u (x N \ t) (12.211) has a similar determinantal structure to the Schur pro¬ 
cess [60]. The Schur process is a probability measure on the sequence of partitions 
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where X (j) := {(A^, • • • , A^jA^ e Z, X? > ■■■ > Xf > 0}, described 
as products of the skew Schur functions s\/^(x i, • • ■ , x n ). For the ascending case (see 
Definition 2.7 in [10]), the probability measure is expressed as 

N N 

Ilrrri- II S\W/\(k-i)(ak) ■ s X (N)(bi, ■ ■ ■ ,b N ), (2.26) 

i,j =l 1 “ ai °i k =l 

where cij, bj, j = 1, ■■■, N are positive variables. We note that s A (fc)/ A (fc-o(afc) is expressed 
as a fcth order determinant and s X (n) (b\, • ■ ■ , Av) as a Nth order determinant by the Jacobi- 
Trudi identity [50j . 

s\hx{xu ’ ’' , x n ) = det (h Xi - N+j -i(x u ■ • • , x n ))J =1 , (2.27) 

where h^{x i, • • ■ ,x n ) is a complete symmetric polynomial with degree k and t'(A) is the 
length of the partition A. Thus (12.211) and (I2.26p have a common structure of N products 
of determinants with increasing size times an TVth order determinant. 

In the following we provide the relations about two marginals of R u (x N \ t ) (12.191) . from 
which (I2.18P immediately follows. For this purpose, we give two formulas for /^(x) and 
/b(x) (I2.20p . First we define a multiple convolution g* ( ' m ’f(x) m — 0,1,2, • • • for a func¬ 
tions f(x) on R and an integral operator g with the kernel g(x — y) as 

/ OO 

dyg(x - y)g* [k ~ 1] f(y), k = 1,2, • • • . (2.28) 

-oo 

Using this definition, the formulas are written as follows: 

Lemma 5. We regard all integrations below as the Cauchy principal values. For (3 > 0, 
a G C with —j3 < Re a < 0 and m — 0,1, 2, ■ ■ ■, we have 

r B ( -> e -=gcot(^))V, (2.29) 

f*B m) fF(x ) = q m {x)f F (x), (2.30) 


where q m (x) is an mth order polynomial with the coefficient of the highest degree being 
1 /ml. 


A proof of this lemma will be given in Appendix [A] The polynomial q rn {x) in (I2.30P is 
defined inductively by (IA.llD - (IA.13p . But in our later discussion we will not use its explicit 
form. 

From (I2.13p and (\2.29b . we readily obtain for m — 0,1, 2 • • •, 


(z; t) = F j+m ,k(x I t), 


(2.31) 


where we define /s(x) := fs(—x). 

Using (I2.30P and (12.311) . we obtain the following relations. 
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Theorem 6. Let the measures dA\ and dA 2 be 


dA 1 = dx!f\ dA 2 = dxP. 


Then we have 


f N 

/ dA 1 R u (x N ]t) = G(xP, ■ ■ ■ ,x[ N) ]t)Y\f F (x { p - u ), 

7 R JV(JV-l)/2 “ 


3 = 1 
N 


dA 2 R u (x N ;t) = W(x[ N \ ■ ■ ■ JJ/f( 

3 = 1 


IN),. 


JR N ( N-l)/2 

Here G(x[ i] , ■ ■ ■ ,x[ N ^]t) is defined by (12.121) and 


(N) s 
Xj — U). 


N—l 


W(x [ N \••• ,4?°;*) = qj - ?j) -det (F y 

j=i 






iV 


j,fc=i 


(2.32) 

(2.33) 

(2.34) 

(2.35) 


where qj(x ) is defined below (I2.30P . 

We easily see that (12.181) can be obtained from these relations (12.331) and (j2.34p : Inte¬ 
grating the both hand sides of them over the remaining degrees of freedom ((44 • • • , 4^' 
for (I2.33P and (4 , ■ ■ ■ ,44) for (I2.34p h we get two different expression about the inte¬ 
grated value of R u (x n ] t ) 


/ R JV(JV+ 1 )/2 


J R N(N+l)/2 


r N 

dx N R u (x N ]t) = / Y\dxPf F (x ( p - u) ■ G{xp. 

Jr n fJi 


i ' > 


{N) -t ) (2.36) 


3 = 
N 


dx N R u {x_ N \t) — / Y[ dx j N) fF( x j N) ~ u ) ' W(x[ N \ ■ ■ ■ ,x^;t). (2.37) 


3 = 1 


where dx N = rii<i<j<Ar X P ■ RHS of the second relation is further rewritten as 

N 


r . i 

/ Yldxpf F (xp-n) • — 

jRN j =1 ' aW&S N 


(2.38) 


( X(jiN) p) , • • • , X a (N) , f) , 

j=i o-WeSiv 

and the symmetrized W(aq, • • • , Xn ; t) in this equation is nothing but W(xp\ • • • , 4^4 t) (12.21) 
since 

1 

m 

a ( N )es N 

„(N) 


E W( - x A N )(l)i - ’ ' j X rr( N ) (V) 
t (n)gs n 


• 4) = m • det (®-i(4' v> )) JJfc , 1 ** (n^rt«)) , l=1 


= ,44^)- 


(2.39) 
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Here in the second equality we used the fact that q 3 (x) is a jth order polynomial with the 
coefficient of the highest degree being 1 /j\ and F\j(x] t) = ipj{x\t). 

Proof of Theorem El First we derive (12.331) . By the dehnition of (I2.19F LHS of (12.331) 
becomes 

N N 

n M®? 5 - u ) • det (fjB k ~ 1)F ij (x[ N ~ k+1) ] *)) . _ • (2.40) 

3 =1 

Here /s(x) is deffiied below (12.311) . Applying (12.311) to this equation we obtain (I2.33j) . 

Next we derive (12.341) . We see that the factor gL 4 2 rii<i<i< 7 v ~ x \ 3 -\ ^) i n Ru( x n] t) (12.191) 
can be decomposed to 

N —1 / N—k N-k-\-l 

dA 2 n f< (*?> - * i> ) = n n ■ n * 

1 <i<j<N k=l \ i =1 j =1 

(2.41) 

and from (|2.30[) the integration of the factor for each k is represented as 

N-k+l 

n dx ? +t - l> n - rr- k) ~ u) 

1 <i<N-k j= 1 

= q N -k(x^l k+ 1 - u)f p (x < £l k+ 1 - u) (2.42) 

where q m (x ) is given in (I2.30p . Eq. (I2.34p follows immediately from this relation. □ 

3 Dynamics of the two marginals 

The purpose of this section is to have a better understanding of the two quantities, 

W{x i,--- ,£v;t) (12. 2p and G(x i,--- ,XN]t) (12. 12p . which arose as partially integrated 
quantities of R u (x N ]t ) (12.19p in Theorem [6] (for W a symmetrization is also necessary, 
see (12.391) ). We will Erst consider the evolution equations of these two quantities. Next 
we will see that the zero-temperature limit of the equation for W(x i, • • • ,Xjy]t) is noth¬ 
ing but the evolution equation for the Brownian particles with reflection interaction while 
W{x i, • • ■ , xn] t ) satisfies the one for the GUE Dyson’s Brownian motion [3T]J regardless of 
the value of (3. Furthermore we will fold that our idea using R u (xn> 0 i n an enlarged space 
]g> n(n+i )/2 ^Theorem [6]) is similar to the argument in [78] although we need a modification 
of [78] about the ordering in an enlarged space. 

3.1 Evolution equations of G{x i, • • • , t) and W(xi, • • • , xy; t ) 

Let us hrst summarize the properties of Fjk(x ; t ) j, k G {1, 2, • • • } (12.131) all of which are 
easily conErmed by simple observations: 

F lk (x;t) = ^k(x;t), (3.1) 
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/ 3 2 


7 T- 


d Id 2 

qI F j k (x]t) = - ^F jk (x;t), 

/‘OO 

/ dxf B (x - y)F jk (x ; t ) = F j+ i fc (j/; t), 

/—OO 

( ^ X j+ 1 e /3(x j+1 -xj) _ Y^j+lk^j+lit) — Fjk (Xj , k ), 


(3.2) 

(3.3) 

(3.4) 


where ipk(x',t) in (13.11) and f B (x) in (13.31) are defined by (12.11) and below (12.31 j) . Eq. (13. 3 1) 
is equivalent to (12.311) while (13. 4 j) is obtained from the relation 


ft 2 r°° e f (y~ x ) 

7r 2 L ^ e^ y ~ x "> - 1 6 



e 


—bx 


(3.5) 


for |Re 6| < j3/2. This relation is easily given by (I2.29jl with a = b - f3/2. 

We see that due to (13.21) and the multilinearity of a determinant, G(x i, • • • , xa q f) (j2.12[) 
satisfies the diffusion equation. 


d_ 

dt 


G(x i, • • • ,Xjv;t) 



(3.6) 


In addition, by (13. 4p . it satisfies the condition 

/3 2 r°° 

“^2 ]_J x ^ eP{xj+1 - Xj) _ ■■■,XN]t) = 0. (3.7) 

for j = 1, 2, • ■ • , N — 1. Though this condition is unusual, we will see that it is regarded as 
a finite temperature generalization of the Neumann boundary conditions at Xj = x J+ \ , j = 
1, • • • , N — 1 in the zero temperature limit (see (13.19ft ). 

On the other hand, from (j3.2[) with the harmonicity of the Vandermonde determinant 
in (12.21) . we see that ,xa r',t) satisfies the Kolmogorov forward equation of the 

GUE Dyson’s Brownian motion [31], which is a dynamical generalization of the GUE, 


d_ 

dt 


w(x 1 , • • • ,XA7 ;t) 



The time evolution equation for the GUE Dyson’s Brownian motion can be transformed 
to the imaginary-time Schrodinger equation with free-Fermionic Hamiltonian (e.g. see 
Chapter 11 in |3B]). On the other hand note that the density function of the Whittaker 
measure (11.9)1 does not solve such a simple free-Fermionic time evolution equation (13.81) . 
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3.2 The zero-temperature limit and a Brownian particle system 
with reflection interactions 


Let us consider the zero temperature limit of the equations (I3.6[) with (j3.7[) and (13. 8p . Note 
that for x ^ 0, 

- lim f B (x) = l>o(z), - lim f B (x) = lim f F (x) = l< 0 (aj), (3.9) 

p—>■ oo p—>-oo p—>-oo 


where f B (x) is dehned below (12.311) and l >0 (:r) and l< 0 (x) are the step functions dehned 
by 


l>o(z) 



x > 0, 
x < 0, 


l<o(®) 


0, x > 0, 

1, x < 0. 


In addition we have 


(3.10) 


lim Fjk(x\ t) =Fj- k {x- 1 t) ) 

/3—i-oo 


(3.11) 


where J r n {x', t) is dehned for n G Z and e > 0 as 




dx e - Xx+xH / 2 

2ni X n 


(3.12) 


Here we summarize a few properties of the function which are the zero temperature limit 
of (I3.13p ~fl3.16p for F jk (x;t). 


e ~ x2 / 2t ( 1 \ i 

XF-k{x\ t ) = lim i) k {x\ t) = - — H k (x/V2t), k = 0,1,2, -- 

0^og y/ 2 nt \ 2 tJ 

d 1 d 2 , . 

ry 

/ dxF n {x\t) = -F n+l (y]t), 


d 

—Fn&t) = -Fn-liyit), 


(3.13) 

(3.14) 

(3.15) 

(3.16) 


where in (13. 13p . ijj k (x]t) is dehned by (12.11) and H k (x) is the kill order Hermite polyno¬ 
mial [5]. The second equality in (13.131) has appeared as (I2.3j) . Note that (13.161) corresponds 
to the zero-temperature limit of (13.4h . since RHS of (13. 5 p goes to —be~ hx in the zero- 
temperature limit and thus the integral operator with the kernel 7r 2 e^^ y ~ x ^ 2 /fd 2 {e^^ y ~ x ^ — 1) 
is equivalent to differentiation in the zero temperature limit when its action is restricted 
to e~ bx . 

Let Q{x i, • • • , Xn] t ) be the zero-temperature limit of G(xi, ■ ■ ■ , Xn] t) (I2.12p dehned on 
M. N . From (13. lip , we hnd 


Q (xi, • ■ ■ i Xn 1 det {Fj~ k (xn—j+i > t)) j k= \ • 


(3.17) 
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The function Q{x i, ■ • • , xj y]t) appeared as a solution to the Schrodinger equation for 
the derivative nonlinear Schrodinger type model im As discussed in ra, using (I3.14p 
and (13.161) with basic properties of a determinant, we find that for x\ ^ ■ ■ ■ ■=£ xn, 
Q(x\, ■ ■ ■ ,XN‘,t) satisfies the diffusion equation, 


d_ 

dt 


Q{x i, • • • ,X N ]t) 



with the boundary condition 


—Q(xi, ■ ■ ■ ,x N ]t)\ 
dxj 


Xj—>Xj + l 


0, for j 


1,'' • ,N — 1. 


(3.18) 


(3.19) 


The probabilistic interpretation of Q(xi, ■ ■ ■ ,xjsr]t) has been given in [75] . Let X t {t), i = 
1, • • • , N be the stochastic processes with iV-components described by 


Xi{t) — yi + Bi(t ) + L i (t), (3.20) 

where i/j 6 I satisfying y x < y 2 < ■ ■ ■ < Vn represent initial positions, Bi(t ) denotes the 
standard Brownian motion and L~(t) is twice the semimartingale local time at zero of 
Xi — Xi _i for i — 2, • • • ,N while L±(t) = 0. The system (13.201) describes the TV-Brownian 
particles system with one-sided reflection interaction, i.e. the ith particle is reflected from 
the i — 1th particle for i = 2, 3, • • • , N. In [78], Warren found that the transition density of 
this system from yi to Xj, i — 1, ■ ■ • , TV is written as Q(x\ — yi, ■ ■ ■ ,xjy — y at; t). Such kind 
of determinantal transition density was first obtained for the totally asymmetric simple 
exclusion process (TASEP) in [7T]. Furthermore, based on the determinantal structures, 
various techniques for discussing the space-time joint distributions for the particle positions 
or current have been developed for TASEP [I51IT7I - (211I561I641I65] and the reflected Brownian 
particle system (I3.20p [321133] ■ 

On the other hand, we have seen in (12.bp that the zero temperature limit of W (xi, • • • , Xjv; t ) (12.2p 
is the GUE density Pgue(xi,-- - ,XN',t ) (11.4p . Note that Pq ue(xi, • • • ,xjsr',t ) also satis¬ 
fies (13.8p since it holds for arbitrary /3 i.e: 


d i N o 2 

ue(xi, • • • , xjv; t) = - ^ ^2-Pgue(xi, • • • , x N : t ) 

3 =1 J 


d 


N 

^ dxj 

3 —1 


( 


N 


\ 


y;— 

m= 1 1 m 

\m^j 


-Pgue(xi, • • • , xn] t). (3.21) 


From (12.60 . (13.9D . and (13. lip , we find that the zero-temperature limit of (I2.18P is 

r N „ N 

/ T\dxfG(x 1 ,---,x N ;t)= / T\dxj ■ Pgue(xi, • • • ,x N ;t). (3.22) 

J (—00,u] N £—1 J ( — 00 % u\ N X_1 


(-OO.u]^ - = 1 
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In [78] Warren showed that this relation, which connects the two different processes, is 
obtained in the following way. First one introduces a process on the N (jV+1) /2-dimensional 
Gelfand-Tsetlin cone whose two marginals describe the above two processes. The Gelfand- 
Tsetlin cone GT^, k — 1, 2, • • • is defined as 

GT fc := {(x (1 \ • • • , x^ fc) )| x® = (x^, • • • ,x^) G M* with i — 1, • • • , k, 

x^ 1 '* < x< x^ +1) with l<£<m<k — 1}. (3.23) 

For the graphical representation of an element of GT&, see Fig. |Tj (b). Next we intro¬ 
duce a following stochastic process on GTjv. Let - ,X^ N \t)) with X^\t) = 

■ ■ ■ ,Xj j \t)) be a process defined by 

+ y® + L^~{t) - L?' )+ (t), 1 <i<j<N, (3.24) 


where B^\t) are the N(N + l)/2 independent Brownian motions starting at the origin, 
yti) re p resen t the initial positions and the process L^' 1 (t) and L^ + (t) are twice the semi¬ 
martingale local time at zero of — X^ ^ and X^ 1 — xjfi( * ^ respectively. Eq. (I3.24p 
describes the interacting particle systems where each xj J> ( t ) is a Brownian motion reflected 
from X^^it) to a negative direction and from X^~ l \t) to a positive direction. In [16], 
Borodin and Ferrari also introduced similar processes on the discrete Gelfand-Tsetlin cone 
where the probability measure at a particular time is described by the Schur process [60]. 

The pdf of the system (13.241) at time t can be given explicitly : For the case of y\^ = 0, 
it is expressed as 


Qgt(x n ] t) 



l<i<7<./V 


n exp 

n- 


k =1 



t k 1 \f 7 2jvt 


Igt(^at)- 


(3.25) 


where x N is dehned above (j2. 19j) and lcT^h-) represents the indicator function on GT fc . 
The pdfs of the two marginals, (x^, • • • , x^) and (x^ , ,xffl) for Q G t(x_n'A) was 
obtained as follows: 


Proposition 7. (Proposition 6 and 8 in (78H J 

dA 1 Q GT (x N ]t) = G(x,x[ N) ]t ) JJ l >0 (x] y+1) -xj j) ), 


N-l 


/ R JV(iV-l)/2 


dA 2 Qgt(x n ] t) = N\P GEE (x[ N \ ■ ■ • ,x{J°;t) JJ 1 >o(x ( j N} - x^J), 


/ R JV(AT- 1 )/2 


3 = 1 


(N).. 


N-l 


, W JN )< 


(3.26) 

(3.27) 


3 = 1 


where G(x i,--- , xn ] t), P G \j E (x { ( W \ ■ ■ ■ ,xffl;t), l>o(x) and dAi , gL4 2 are defined by (13. 1 71) . 
(11.41) . (13.1 Op and (12.32P respectively. 
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Remark. Note that Q{x^\--- , xf > ]t) in (13,26j) can be replaced by an arbitrary func¬ 
tion on M ,v such that it corresponds to ,x^;t) in the region xj'^ < x^ 21 < 

• • • < x^. For later discussion on a generalization of finite temperature, we chose it as 
Q{x^\ • • • , xf^; t) on the whole R N . 

We see that the relation (13.221) is obtained from this theorem. By decomposing the 
integral on x N in two different ways, we clearly have 


N 


'(-oo.up^- 1 )/ 2 


dx N Q GT (x N ]t) = 


'(-00 ,u] N - =1 


J^J dx{ 


U) 


l R N(N-l)/2 


dAiQa^x^t) 


N 


U dx 


(N) 


J(-oo,u] N j =1 

Applying (I3.26|) and (13. 2 7ft to this equation, we get 
„ N N -1 


l R N(N-l)/2 


dA 2 QGT(x N ]t). (3.28) 


'(-oo,u] N J =1 


Y[dx[ j) g(x[ 1] , ■ ■ • ,4 jV) ;f) JJ 1 >0 (xS j+1) -xj j) ) 


3 = 1 


,. N N -1 

/ II dx< j N) N\P G ve(x[ N \ • • • , Xfiprf) JJ l>o(x 

J(-oo,u] N f=i f=l 


(TV) _ (N) % 

X 3+A 


(3.29) 


Due to the symmetry of Pgue(^i,"‘ ,%n'A) under the permutations of Xi,--- ,xjv, we 
readily see that RHS of this equation is equal to RHS of (13.221) . Also we find that LHS 
of (13.29(1 becomes 

JV-l 

.b+i) JJh 


[ dx S j) ■ £(xi 1} , • • • , xS t) JJ 1> 0 (4 } - x { { 

j=1 j= i 

„ N JV-l 

= / n dx[ j) ■ Q{x S 1} , • • • , x[ N) ]t) n (l>o(ajS J+1) - } ) + 1 >o(^S j) - +1) ) 

J{-oo,u] N t=l 7=1 V 


J=1 

N 


Yl dx i ] •£(x! 1) ,--- ,x^;t), 


R-oo.u]^ J =1 

where in the first equality we used for k — 2, ■ • • , AC and (xf \ • • • , x^) G (— oo, u] N 


(3.30) 


k—2 


dx[ k) g(x[ L \ ■ ■ ■ ,x[ N) ]t) IJl>o(xi : 

'(-oo,«l J =1 


0+1) Rib 


x ">) ■ l >0 (x[ k ~ 1] - xf } ) = 0. (3.31) 


Note that Q{x^\ • • • , x^; t) is defined on R N and is finite even outside the region x^ < 
xf^ < ••• < x[ N) . (See Remark, of Proposition [TJ ) Eq. (13.311) is obtained from the 
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following observation: putting the last factor l> 0 (4 fc 1 ' ) — xf^) in the N — k — 1th row of 
the determinant , • • • , x^; t ) in (13.311) then applying (I3.15p . we get the determinant 

which has the same two rows. 

Thus (I3.22p is obtained from Proposition [3 This is similar to the situation of (I2.18P and 
Theorem [6j This naive observation gives us the impression that the pdf QgtGLaO t) (13.25P 
is the zero-temperature limit of the weight R u (x N : t) (12. 19p . However in fact this is not 
the case. Let Tl u {x_ N ]t) := lim^oo R u (x N ] t). From (13. 9p and (13.111) one has 

R-u{x N \t) = (-l) JV(JV - 1 ) /2 d e t (j r l - i (x ( f ) -,t)^ ,_ i l>o( x f-i ) ~ x j k) )- (3-32) 


From (12.5p and (13.131) . it is further rewritten as 


Ru(x N ]t) 



1 v N (x N ), 


(3.33) 


where 1 v k (x k ) i s ^ ie indicator function on an ordered set 14 defined by 
14 := {(x (1) , • • • , x^)| x^ = (x[ j \ 


,x^) e < x*j 


(m+l) 




1 < 




< m < k — 1}. 

(3.34) 


For the graphical representation of an element of (I3.34p . see Fig. CD (a). Comparing (I3.25P 
with (I3.33p . we see that they have the same form but their supports (GT^r and Vn) are 
different. We further notice that Vn with an additional order x^ < xf n+1 ^, 1 < l < m < 
N — 1 corresponds to GT]y. 

Hence our approach using TZ u (x N ]t) can be regarded as a modification of Warren’s 
arguments on GTat to the ones on the partially ordered spece Vn- Let us focus on two 
marginals (x^, x[ 2 \ ■ ■ ■ , x^) and x^\ ■ • • ,x^) for 7 Z u (x N ;t) (13.331) . By taking 
the zero-temperature limit of Theorem [HI we have the following analogue of Proposition [3 


Proposition 8. 


/ r jv(;v— 0/2 


N 


dAi lZ u (x N ] t) = ,x[ N) ;t) JJl >0 (n -x^), 


(N).. 




(3.35) 


3= 1 


r N 

/ dA 2 TZ u {x N ]t) = P u (x^\--- ,x^ ) ;t) TTl>o (u-x ( j N) '\ , (3.36) 

J r jv(jv-i)/2 v j=i V 

where for the definition of dA\ and dA 2 , see (12.321) . Qfx[ \ ■ ■ ■ ,x^;t) is given by (13.171) 
and 


=n o-i)J- • n 

j =1 ’’ l<j<k<N 


N -1 x (N) 


JN) _ (AT) 

Xj x k 


n 

3 = 1 


\ 2 

) / 2 1 


18 


y/2nt 

(3.37) 




































Proof. It is obtained by taking the zero-temperature limit ({3 —> oo) in Theorem [6] □ 

As discussed in (I2.39p . Pgue(^i 5 • • • ,%N]t) (II .4ji can be interpreted as the symmetric 
version of P u (x i, • • • , x n] t ): 

Y (lb'" , x <r(N)( N y,t) = Pg\je(x[ N) , ■ ■ ■ ,4?°;f). (3.38) 

■ crWeSN 

Therefore by the similar discussion in (I3.28p . we see that the relation (I3.22p is obtained 
also from Proposition [HJ 

The fact that both Proposition [7] and [H] lead to (I3.22p implies the relation 


/ R JV(AT+ 1 )/2 


dx N lZ u {x N \t) = 


l(-oo,u] N(N+1)/2 


dx N Q GT (x N ]t). 


(3.39) 


This equivalence of their integration values is generalized in the following way. 


Proposition 9. Let / sy m(2iv) the function defined above fl2.22p . Then we have 


J R N(N+l)/2 


dx N f^ ym (x_N )‘P'U (d-N * f ) 


'(-oo , u ]JV(JV+l)/2 


dx N f sym (x N ) Qgt(x n ] t ) (3.40) 


An essential step of the proof of this proposition is represented as the following 

Lemma 10. 

Y sgna (N) l VN (xfi N ) = Y sgnu (Ar) l G T(aZjv) (3-41) 

aU)eSj,j=l,- ,N aU)eSj,j= l,-" ,N 

The proof of this lemma will be given in Appendix [Bl Using this lemma we readily 
derive Proposition [HI 

Proof of Proposition [9j Substituting the definition of TZ u (x N ]t) (13.331) into (I3.40p . we 
see that the LHS of (I3.40p is rewritten as 


/ R iV(JV+l)/2 


dx N f sym (x N ) U 1 >o(u-4 /l, ) e ^ /2t ■ n 


k =1 


^Ar/symfejv)n i> 0 ^ _ 2 : l fc) ) e ^ ^ ^ ' II 

k=l 1<2<7 < 


X W _ X W 


l<2<_7<Ar 

x Y sgna (Ar) ly(^) 




/ R JV(AT+l)/2 


X (N) _ X (N) 


l<i<j<N 

x Y sgna (Ar) l G T(^) 


N 


aU)eSj,j=X,-,N 

dx N f sym (x N )t\l >0 (u ~ X { 2) Y QGTixfif-fi) 

• /RjV(Af+1)/2 *=1 aU)e S] , 3 =l ,-.,N 

where in the second equality we use Lemma ITOl 


(3.42) 

□ 
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4 Fredholm determinant formulas 


4.1 A Fredholm determinant with a biorthogonal kernel 

The function W{x 1 , • ■ ■ , apv;t) (11.611 has a notable determinantal structure that it is de¬ 
scribed by a product of two determinants. This allows us to apply the results of random 
matrix theory and determinantal point processes developed in | 43ll76j and to get the Fred¬ 
holm determinant representation. 

To see this we provide a lemma. Let (f)j(x ; t ), j = 0,1,2,--- be 


(j)j(x]t) 


1 

2ni 


dv e vx ~ v2t/2 


T(l + v/P) N 
lb+ 1 


(4.1) 


where the contour encloses the origin anticlockwise with radius smaller than (3. We find 
<j>j{x-,t) and ^ k (x]t) (12.11) satisfy the biorthonormal relation: 

Lemma 11. For j, k G (0,1,2,---}, we have 


dx 4>j(x] t)ip k (x-1) = 8 j}k . 


(4.2) 


Proof. Substituting the definitions (12.ip and 

poo 

/ dx 4>j(x] t)ijj k (x; t) 


into LHS of (14.21) . one has 


(27r) 2 i 


dv 


dw e 


-i„wit/2 f r (i +v/p) \ ri (»«)* r dxe ^. 


T(l + iw//3)J vi +1 J_ 


(4.3) 


As the integrand in this equation is analytic on C with respect to w, we can shift the 
integration path as w = w' — iv, w' G M. Then using 


1 

27T 



dxe {v ~ iw)x 


1 f°° 

— J dxe- iw ' x = 8{w') 1 


(4.4) 


we find 


dx (j)j(x] t)if) k (x] t ) 


1 

27 t i 


dvv k j 1 = 5j }k . 


(4.5) 


□ 

The residue calculus shows that the function (f>j(x ; t) is a jth order polynomial in x and 
the coefficient of the highest order is 1/jl. As the Vandermonde determinant in (12 .2 [) is 
expressed as 

JJ (x k - Xj) = det (4 _1 )^ fc=1 = det (t7 “ , ( 4 -6) 
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(4.7) 


W(x i, • • • , Xn ; t) is rewritten as a product of two determinants 


W(x i, " ' > x N’,t) = det i(x fc ; t))^ =1 det (V’j-ifc*; t))£ fc=1 . 

From Lemma EH and (14.71) . we obtain a Fredholm determinant representation for the mo¬ 
ment generating function. Throughout this paper, we follow [TDj for the notation on 
Fredholm determinants. 


\N 


Proposition 12. 


N 


]^[ dxj g(xj) ■ W Oi, • • • ,x N ;t) = det (1 - gK) 


L 2 ( R) 


(4.8) 


' j =i 


where g{x) is an arbitrary function such that the left hand side is well-defined and in the 
right hand side det (1 — gK) L2 ^ represents a Fredholm determinant defined by 

det (1 - gK) L2{m = ^ —— / J ~[dx j g(x j ) • det {K(x h x m ]t )) k lm=1 . (4.9) 

k =o ' j =i 

Here g{x) = 1 — g(x) and K(x,y;t ) is written in terms of the biorthogonal functions 
ipj(x,t) (12.11) and <fk(x,t) (14.11) as 


N -1 


K(x,y,t) = Y. i j> k (x\ t)ip h (y]t). 


(4.10) 


k =0 


Proof. We readily obtain this representation by applying the techniques in |76j with 
Lemma EH to LHS of (14. HD - For reference, here is an outline of the proof. First, using the 
Andreief (Cauchy-Binet) identity (12.161) . we have 


N 


Y[d Xj g( Xj ) -W(x i,-- - ,x N \t) = det ( j^dxg{x)(f) j - 1 {x\t)'if k _ 1 {x- 1 t) 

3 =i 


N 


j,k= 1 


= det I / dxfa-ifattyk-i&t) - / dxg(x)<pj-i(x;t)i(3k-i(x;t) ) 

\«/R J M / 

= det (<5j- fe . — A jt k)f tk=1 , (4.11) 

where Aj jk , j, k — 1, • ■ • , IV is defined as 

A** = / <fo<K®)0j-i( a; ;*)V’fc-i(®;f)- ( 4 - 12 ) 

In the first equality of (14.111) . we used (14.71) with (12. 16jl and in the last one we used 
Lemma EH We further rewrite A ]k as 


Ajk = / dx B(j, x)C(x, k) 


(4.13) 
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by using 


B(j, x ) = <t>j_ i(s; t), C(a;, fc) = £(®)-0jfc-i(®; t). (4.14) 

Applying the identity for Fredholm determinants, 

det (8 j)k - A jtk )* k=1 = det(l - BC) L 2 ({1j2j ... !N}) = det(l - CB) l2{wl) , (4.15) 

and noting 

JV-1 

( CB)(x,y) =g(x)^2<i>k{x-,t)il> k (y,t), (4.16) 

fc =0 

we arrive at our desired expression. □ 

Combining this proposition with Theorem [21 we readily obtain 

Corollary 13. 

/ e-^Z N (t) \ 

E ( \ = det (! _ f u K) L , m (4,17) 


where the right hand side is the Fredholm determinant (j4.9[) with the kernel f u (xi)K(x Xj\t), 
fuixj) = 1 — fF^Xj — u), and K(xi,Xj]t ) is defined in (j4.10j) . 

As in (12.311 . we see 


lim (j) k (x;t) = -— : ® dv 


p— KX> 


2ni 


0 vx—v 2 t/2 
n lfc+1 


k 

= ^l 


k\\ 2 


3), 


Vzt) 


(4.18) 


which is due to another representation of the nth order Hcrmite polynomial H n (x) (see 
e.g. Section 6.1 in [5]), 


H n (x) = 


n\ 

2ni 


dz ■ 


Zlxz—z 1 


?n+1 


(4.19) 


where the contour encloses the origin anticlockwise. From (12. 3 I) and (14.181) . we hnd 


lim K(xi,X 2 ',t) 

( 3^-00 


e -A/ 2 t (VI H k { Xl /y/2i)H k {x2/\/2i) 

V2^i ^ 2 k k\ 


(4.20) 


Here RHS appears as a correlation kernel of the eigenvalues in the GUE random matri¬ 
ces [52] . 

Thus K(xi,Xj]t) is a simple biorthogonal deformation of the kernel with Hermite poly¬ 
nomials which appears in the eigenvalue correlations of N x N GUE random matrices. 
Using this Fredholm determinant expression (14.1 7H . we can understand a few asymptotic 
properties of the partition function by applying saddle point analyses to the kernel as will 
be discussed in Sec. [5j 
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4.2 A representation from the Macdonald processes 

In [57], O’Connell first introduced the probability measure on M. N which is called the 
Whittaker measure m(x i, • • • , xn', t) fljLi dxj whose density function m(x i, • • • , x^; t ) is 
defined in terms of the Whittaker function v I'a(^i, • • • ,xn) (see [57]), 

m t (x i, ,x N ;t) = ^o(/3xi, • • • , f3x N ) / dX i, • • • , /3x N )e^ =1 A ^ /2 s N (X//3), 

J(iR) N 

(4.21) 

where throughout this paper we denote A = (Ai, • • • , Ajv) and sn( A) is defined by (12.101) . 
Then he showed the following relation about the distribution of the free energy F N {t) = 
— log(Z N (t))/(3 (see Theorem 3.1 and Corollary 4.1 in [57]). 

Prob (—F N (t) + — ^ log/3 2 < s\ = f dx i f Y\ dxj ■ 7n(xi, ■ ■ ■ , x N ;t). (4.22) 

VP/ /(-oo,s] J l^- 1 j=2 

The density function m(x i, • • • ,XN]t) (14.211) is also a finite temperature extension of 
-PgueO^Ij •'' ,XN',t) (11.41) . Actually it has been known that m(x i, • • • , xjv;f) converges 
to T’gue(^i, • • •, xn] t) in the zero-temperature limit. (See Sec.6 in [57]). In contrast to 
W[x i,--- ,XN',t) (12. 2p . however, this extension does not inherit the determinantal struc¬ 
ture which Pgue(^i, ■ • • ,%i v',t) has and thus we cannot apply the techniques in random 
matrix theory which is useful especially for asymptotic analyses of the GUE. This fact 
necessitated the developments of new methods p fIUl[TTlU^l[m[211l291l3Ul[55IIB^] . By using 
the techniques of the Macdonald difference operators m and the duality [H], one can get 
a Fredholm determinant expression for the moment generating function of the partition 
function, which allows us to access the asymptotic properties. 

Proposition 14. ( [TO/) 

E e l»rir = det (1 + L) 1I(Co) (4,23) 


where C 0 denotes the contour enclosing only the origin positively with radius r < (3/2 and 
the kernel L(v,v'-t) is written as 


1 f 7 r / w N e wH/ 2 -wu x r(l + v'/(3) N 

2m J iM _ +s sin(n' — w)/(3 v ,N e v ' 2t l 2 ~ v ' u w — v T(1 + w/(3) N 
Here 5 satisfies the condition r < S < (3 — r. 

We can show the equivalence between the two expressions (14.171) and (14.231) . 

Proposition 15. 


(4.24) 


det(l — fuK ) X/ 2 (m) — det(l + T)l 2 (c 0 ) (4-25) 

where f u (x) = 1 — [f(x — u) and K(x,x';t) and L(v,v';t) are defined (14.101) and (j4.24[) 
respectively. 
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Proof. 

Substituting the definitions (12.ip and (14.ip into (14.101) . we have 


K(x, x'; t ) = 


N—l 


dv i 

Co J zR+<5 


dw e 1 


z- wx '-( v i- w i)t/2 r (! + v/P) N 1 ^ 


T(1 + w//3) N v 




k =0 


= ® dv 


dw e 1 


Co J zR+<5 


c-wx'-(v 2 -w 2 )t/2 r(l +V//3) 1 - (w/v) 

T(1 + w//3) N v — w 


\N 


(4.26) 


For the definition of Co, see below (14.231) . Here we changed w —> —iw in (12. ip and shift 
the path of iv by 5 which is larger than the radius of v. We notice that although the last 
expression in (I4.26j) consists of two terms proportional to l/(v — w ) and (■ w/v) N /(v — w), 
the integration of the term proportional to 1 /(v — w) with respect to v vanishes. Thus we 
see 


f u {x)K(x, x') = -f u (x) 


dv / 

Co J iR+<$ 


dw 


e vx-wx'-(v 2 -w 2 )t/2 / pQ _|_ v jp\ 


W — V 


W 


T(1 +w//3) V 


N 


= — ® dv A{x, v)B(v, x') 
Jco 


(4.27) 


where we set 


A(x,v) = (hi 


B(v, x') = 


dw 


- > —wx'+w 2 t/2 


N 


W 


w — v \T(l + w/f3) 


N 


(4.28) 

(4.29) 


Here we use the relation for Fredholm determinants, det(l — AB) L 2 ^ = det(l — BA) L 2 ^ Co ^, 
where the kernel — {BA){y,v') on RHS reads 


dx B(v, x)A(x, v') 


dw 


e (w 2 —v' 2 )t/2 / W Y{1 + v'/j3) 


w — v \dT(1 + w//3) 


N 


dxf u (x)e^ v '- w)x . 


Using the relation 


dx ■ 


7r 


1 + e x sin 7ra 
we perform the integration over x in (14.301) as 

/ * r>n _gf3(x-u)-\-(v'—w)x 


, for 0 < Re a < 1, 


dxf u [x)e 


(v'—w)x _ 


0 (v r —w)u 


dx ■ 




1 + eP( x A sin \{y' — w)tt//3\ 


(4.30) 


(4.31) 


(4.32) 
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Note that because of the conditions 0 < r < (3/2 and r < 8 < (3 — r (see below (I4.23P 
and (14.24p respectively), (14.31 p is applicable to the above equation. Thus from (I4.30P 
and il 1.321) . we have 


dx B(v, x)A(x, v') 


1 f ir/P w N e wH/2-wu 1 r{l + v'//3) N 

2 m J is _ + s sin [(?/ — w)n//3] v' N e v ' 2t / 2 ~ v ' u w — v T(1 + w//3) N 
L(v,v';t) (4.33) 


□ 


5 The scaling limit to the KPZ equation 

In this section, we discuss a scaling limit of the O’Connell- Yor polymer model. When both 
N and t are large with its ratio N/t fixed, it has been known that the polymer free energy 
Fjsr(t ) defined below (II.ip is proportional to N on average and the fluctuation around the 
average is of order N 1 ^ 3 [531 [72]. Furthermore recently it has been shown in [TO] that 
the limiting distribution of the free energy fluctuation under the N F 3 scaling is the GUE 
Tracy-Widom distribution m This type of the limit theorem has been obtained also for 
other models related to the O’Connell-Yor model |6[[l3l[26[[3ll[58ll77] . These results reflect 
the strong universality known as the KPZ universality class. 

Although we expect that the same result on the Tracy-Widom asymptotics can be ob¬ 
tained from our representation (14.17p . we consider another scaling limit where the partition 
function goes to the solution to the stochastic heat equation (SHE) (or equivalently, the 
free energy goes to the solution to the Kardar-Parisi-Zhang (KPZ) equation). This scaling 
limit to the KPZ equation has also been known to be universal although in a weaker sense 
compared with the KPZ universality stated above [T1I91I27]. The height distribution of 
the KPZ equation has been obtained for a droplet initial data in [2 U66H69] . Since then, 
explicit forms of the height distribution have been given for the KPZ equation and related 
models for a few initial data [TUfH21123113511391H§11521153] . In particular for the O’Connell- 
Yor model (II.ip . the limiting distribution of the polymer free energy has been obtained by 
applying the saddle point method to the kernel (14.24p [TOUH] . 

In this section, we confirm that a similar saddle point analysis can be applicable to our 
biorthogonal kernel (14.10p . Since our kernel has a simple form, we find that the nontrivial 
part of this problem reduces only to the asymptotic analyses of the functions ipk{%\ t) ( 12 . 11 ) 
and <pk(x;t) ([Tip . 

5.1 The O’Connell-Yor polymer model and the KPZ equation 

Before discussing the saddle point analysis, let us briefly review the scaling limit to the 
KPZ equation. Hereafter we will write out explicitly the dependence on /3 of the polymer 
partition function ( 11 . 11 ) as Z N ^(t). 
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Let Zj t p(t) := e t ^ 2 Zj t p(t), j = 1, • • • ,N. By Ito’s formula, we easily find that it 
satisfies the stochastic differential equations 

dZ ji/3 (t) = Zj_ lt p(t ) - Z^tfj dt + /3Z jt p(t)dBj(t), (5.1) 

where we set Z 0j p(t) = 0 and interpret the second term of this equation as fto type. Now 
let us take the diffusion scaling for (15. 1 [1 : we set 

t = TM, N — TM — XVM (5.2) 

and at the same time we scale fd as 


jd = M~ 1/4 , (5.3) 

then take the large M limit, the scaling exponent —1/4 in (15.3ft is known to be universal: it 
characterize the disorder regime referred to as the intermediate disorder regime [T], which 
lies between weak and strong disorder regimes in directed polymer models in random media 
in 1 + 1 dimension. 

This scaling can be explained in the following heuristic way. Let Bj(t), j = 

1, • • • ,1V be N independent one dimensional standard Brownian motions. For N \, N 2 G 
{1, 2, • • • , N}, we have 

(B Nl (t)B N2 (t)) = tSN 1: N 2 , (5.4) 

where (•) represents the expectation value with respect to the Brownian motions. Now we 
consider its large M limit under the same scaling as (15.2ft i.e. t — MT and 

N k = TM - X k yfM, k = 1, 2. (5.5) 

Noting that liniM^-oo VM5n u n 2 = — X 2 ) under fj5.5f) . we see 

lim M” 1,/2 (B Nl (t)B N2 (t)) = T5{X l - X 2 ). (5.6) 

M-> oo 

This suggests in a heuristic sense, 

lim M~ 1/A B N {t)= [ dsr)(s,X k ), A; = 1,2. (5.7) 

M—>-oo J Q 

Here r/(T,X) with T > 0 and X 6 1 is the space-time white noise with mean 0 and 
5-function covariance, 

(V(T,X)) = 0, ( V (T,X) V (T',X')) = 5(T - T')5{X - X'). (5.8) 

Thus considering (15. 7p . we choose the scaling of (3 (15.31) . 
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Under the scaling (15.2ft and (15.311 . the following limiting property is established. 


lirn VmZnM) = Z(T,X). (5.9) 

M—>oo ’ 

Here Z(T,X) is the solution to the SHE with the 5-function initial condition, 

§fZ{T, X) = 1 --^Z(T,X)+ri(T,X)Z(T,X) (5.10) 

Z(0,X)=S(X), (5.11) 

where rj(T,X ) is the space-time white noise with mean 0 and 5-function covariance (15.8p . 
The SHE (15.101) is known to be well-defined if we interpret the multiplicative noise term as 
Ito-type [81155]. Using this equation, the solution to the KPZ equation can be defined via 

h{T,X)=log(Z{T,X)), (5.12) 


which is called the Cole-Hopf solution to the KPZ equation. Recently a new regularization 
for the KPZ equation was developed in [37] (see also [38]). 

According to ra, a rigorous estimate about the convergence to the SHE (15. 9 p has 
been obtained for the O’Connell-Yor model [53] based on the results in |TJ. This type of 
convergence has been discussed also in interacting particle processes PEI]. For reference 
we offer a sketch of the derivation of (15.91) . For this purpose, we provide the following 
lemma, 

Lemma 16. For Z N ^{t) defined above (15. ip . one has 


OO 

z«At) = E 3 * 

k =0 


l<Ni<-<N k <N ' 


1 h ,i; 

^(o,t) j=1 


N. 


(*;) 


k+1 

n p °^> 

3 = 1 


tj~i 


Nj 


Nj- 1 ) (5.13) 


where Po(t,n) := e~H n /n\ denotes the Poissonian density and N 0 = l,Nk+i = N,sq = 
to = 0, sn = tk +1 = t. A n (s,t) denotes the region of the integration s < t\ < ■ ■ ■ < t n < t 
and the ltd integrals on RHS, referred to as the multiple ltd integrals are performed 

in time order (i.e. the order oft 1 , • • • ,Uvj. 


Proof. By the definition of Z^{t) (II.ip . we have 

N -1 N 

z N ,p(t)=e~ t n^-n« 

J 0<Sl<---<SN-l<t j — l j — l 


1)-- 


(5.14) 


with sq= 0 and the integrand of RHS is expressed as 


N i)\ N 


3 = 1 




np 

3 =1 


P\ B i( s i)-Bj(sj- 1 )-- 
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e e n 

m= 0 l<Mi<---<M m <N j= 1 


. /3( s M 

P( BMjiSMjl-BMjtSMj-l) - - 2- - - 


- 1 


(5.15) 


Here we use the relation on a one-dimensional standard Brownian motion B(t)\ one has 
for t > s > 0 and (3 > 0, 


= £ pn f TT dB(fjh 
n =0 >/A n (s,i) J=1 


(5.16) 


where the Ito integrals on RHS, referred to us the multiple Ito integrals, are performed in 
time order (i.e. the order of f 1; • ■ • ,t N ) [40(51]. Using this, we get 


N P(BjM-Bj 0{aj ^ ^ ~ 


n 

3 = 1 


m oo 


E E IIE*’ 


oo oo 


E^E E E n/ 

fc=0 m=0 rei,-,n m =l j=l J An j ( s m 3 -i,sm 3 ) ^ =1 

filH—+n m =fc 


=0 j=l nj = l J £=i 

nj 

n M i 


m 

oo m 


|^J jt e) 

(5.17) 


Substituting this into (15.141) . and performing the integration on si, • • • , Sjv-i, we have 

m n j 

n o dB M j 


oo oo 


oo p m ,l 3 

z N , P {t) = [ _ 

niH- \-n m =k 


x e 


-t 


m+l 

n 

3 =1 


(•'/,' •+ i)! 


(5.18) 


where we set M 0 = 1, M m+ 1 = IV. Now we introduce the new variables Nj, tj, j — 1, • • • , k 
by the relation 

N ni -^ \-rij—\+l , tnH fn.y_i tMj,l for L 1, • • • , Tlj , J 1, , Ui. (5.19) 

Then one has dB Mj {t Mj ,i ) = d5jv ni+ ... + „._ 1+f (in 1 +-+n ;# _ 1 +^) leading to 


nn-«w* 

i=ir=i 


Mil) 


Y\_dB N .(tj). 

1=1 


Further from (15. 19f) . we have 


—t 


m+l 

n 

3 = 1 


(+ + \Mj-Mi_ i fe+1 /. . \Nj—N,-\ 

_ tt- tj-v 1 1 


(Mj - M 3 -_i)! 


=n« 

1=1 


(JVj - JVj-i)! 


(5.20) 


(5.21) 












where we set N 0 = 1, N^+i = N. Substituting these (I5.20p and (I5.2ip into (I5.18P and 

■voo 

jm =0 ^-^l<Mi<---<Mm<N ■ 


noting the summations <■■■<¥ <n Yl°rn,--- ,n m =i can be summarized as the 


simple form ^2i< Nl <...< Nk < N , we obtain (15. 13j) . 


niH- \-n m =k 


□ 


Note that under the scaling (15 .20 . the Poissonian density Po(t, N ) goes to the Gaussian 
density g(T,X) = exp(— X 2 /‘IT') / \/2nT, i.e. 


lim VMPo(t , N- 1) = g(T, X). 

M-> oo 


(5.22) 


Furthermore by Theorems 4.3 and 4.5 in [f], for a function f(ti, ■ ■ ■ ,tk, Ni, ■ ■ ■ , N^) that 
converges to f(iti, • • • , u^, yi, ■ ■ ■ , y^) under the scaling t t = tqM and Ni = UiM—y^y/M^ i = 
1, • • • , k, we have 


lim „ ro , ,, 

M —>oo j\/f3fc/4 


E 


dBw^tj) ■ f(t\, ■ ■ ■ ,tk‘, Ni, • • • , Nk ) 

i=i 

/ TT du j" / TT d %' ' TT ,Uk‘,yi,--' ,Vk) (5.23) 

^(0;T) = ! ' Jr^ = 1 m=1 


l<Wi<--.<V fc <Af J A k(0,t) j =1 
k n k k 


where rj(t, y) is the space-time white noise with the 5-covariances 115. 81) . Thus from (15. 13p . f!5.22p 
and (15.23H . we have under the scaling (15. 2p . 

lim VMZ NtP (t) 


M—>oo 


lim VhdM 1 / 4 )* 


k =0 


M 3fc / 4 


e /, n^.fe) 

l<JVi<-<JVj.<iV J A k(0,t) j =1 

k +1 

x JJ M 1/2 Po(tj - A'j - .V,- ,) 


i=i 
fe+i 


n n. n rv n. I J- 

e _ n dtj ■ n dy j n ^ ■ n ^~ 

fc=0 ^ A fc( T ) j=l “' Rfc j=l m=l £=1 


-i, Ue. — ye- i), 


(5.24) 


where t 0 = 0, tfc+i = T, y 0 — 0, ?/fc + i = X. Since we easily find that RHS of this equation is 
the solution of the SHE with 5-function initial data (15. 10p . we obtain (15.9p . 


5.2 The asymptotics of the kernel 

In [TO], Borodin and Corwin discussed the asymptotics of the Fredholm determinant (I4.23P 
under the scaling limit to the KPZ equation, especially the limiting property of the ker¬ 
nel (14.241) based on the saddle point method. Here we check that a similar saddle point 
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method is applicable to our biorthogonal kernel (14.101) . The scaling limit we consider is 
(15.9p discussed above, but here we adopt its rephrased version stated in [10], 


.. z N p = i(t = Vtn + x ) 

hm ———- 

N^oc C(N , T, X) 


Z{T,X), 


where C(N, T, X) is 


C(N, T, X) := exp 



Vtn+ x 
2 




(5.25) 


(5.26) 


which is more suitable for our purpose. To see the equivalence between (15.91) and (|5.25j) . 
we rewrite the relation (15. 9 p as 


hm B 

N—too 


-2 


Z N , p {t) = Z(T,X), 


where we scale t, j3 as 


t — N + X 





(5.27) 


(5.28) 


Furthermore focusing on the scaling property of the partition function Z N ^{t) — ZN,i(P 2 t)//3 2 ^ N 1 , 
we hnd 

P~ 2 Z N ,p(t) = p 2Net+ p2 t/2 Z N ,i {P 2 t) (5.29) 

in distribution. Noticing under the scaling (I5.28P 

/ 3 2 t = Vtn + x, p 2 AT e t+ p 2 t / 2 = c(n, t, x), ( 5 . 30 ) 


where C(N,T,X ) is defined in (15 .2fj[) . we hnd that (I5.27P is equivalent to (15.251) . 

For the moment generating function, (I5.25P implies 

hm E (^-^ZnVVtn+x)^ _ E V-e-'z(T,x)^ _ E ^ (5.3i) 

where on LHS, u is set to be 

u — u’ + log C(N, T, X), (5.32) 

with C(T,N,X) (I5.26P . and in the last equality in (I5.3ip we used (15. 121) . The notions 
of the KPZ universality class tell us that the fluctuation of the height h(T,X) and the 
position X are scaled as T 1//3 and T 2 / 3 respectively for large T. Considering them, we set 

h (T, 2 7 }Y) = ~2| + 7 t(HT, Y) - i' 2 ), (5.33) 
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where 7 t = (T/ 2) 1//3 . The first term — 7 j ’/12 = —Tj 24 represents the macroscopic growth 
with a constant velocity. The height fluctuation is expressed as h(T, Y ) and the term Y 2 
reflects the fact that the SHE with delta-function initial data in (15.111) corresponds to the 
parabolic growth in the KPZ equation [211661169] . Thus substituting u' = 7 tS— 7 ^/ 12 — 7 tH 2 , 
X = 27 jY into (15.321) . we arrive at the modified scaling 


7 ^ Ar VTN + 2'ylY „ UV N , T 

« = 7ts - — - 7t>' + iV +---+ 2 7t 1 Y + y lo S y • 


Hence (15.31 j) is rewritten as 


lim E 

TV—>00 


^ e -e—Z NA (VTN+2Y;Y)^ = £ ^_ e 7t(i(T,r)-.))^ . 


(5.34) 


(5.35) 


with the scaling (15.341) . This is the scaling limit of the moment generating function from 
the O’Connell-Yor polymer to the KPZ equation. 

It has been known that RHS of this equation can be represented as the Fredholm 
determinant [24 |1291I30 ], 

E (e- e7T(MT,i ° s) ) = det (l - /C K pz) i2(K) , (5.36) 

where the kernel /Ckpz(£i>£ 2 ) is expressed as 

p 7t(?i-s) r°° 

/c kpz (6,6) = e7T ( gl _ g ) + 1 j d\ Ai(£i + A)Ai (£ 2 + A). (5.37) 

Note that Y does not appear in RHS of this equation. This kernel first appeared in the 
studies of the KPZ equation for the narrow wedge initial condition [21 1661 - 169] . From the 
relation (15.36ft we readily get the distribution of the scaled height h(T,Y) given in (15.33ft . 

By combining the formula (14.171) for the O’Connell-Yor polymer and the limiting rela¬ 
tion (I5.35P from the O’Connell-Yor polymer to the KPZ equation, we can obtain (15.361) by 
showing 


Jdydet 0 - - fn K ) L 2 (R) = det (1 - /C K pz) l2(r) 


(5.38) 


under (j5.34[) . This was indeed already discussed in [TO] by using the kernel (14.241) . Here we 
show that the kernel (15.371) appears rather easily from the scaling limit of our biorthogonal 
kernel ( 14.171) . Using the saddle point method, we get the following: 


Proposition 17. 

lim f u ( Xl )K( Xl ,x 2 -VTN + 2^ T Y) = e^-^XKPz^^). (5.39) 


Here the kernel is expressed in terms of 4>k(x\\t) and il>k( x 2 ',t) defined by (14.ip and (12.ip 
respectively as 

gn —u yy 

f u (x 1 )K(xi,x 2 ;t) = xi _ u </>k(x 1 ; t)Mx 2 ; t ), (5.40) 

+ 1 k =0 
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and we set u to be (I5.34p and 


t 7 T 

Xi = It ii - ^2 


7t F 2 + N + 


(T7V) 1 / 2 + 2 7 2 U 


9 N N T , 

+ 2 7t 1 y—+ yl°g—. (5.41) 


Since the factor in (I5.39j) does not contribute to the Fredholm determinant, we 

get (15.38(1 (though for a complete proof one has to prove the convergence of the Fredholm 
determinant itself, not only the kernel). Note that (15 .40(1 has a similar structure to the 
kernel (14.201) in the GUE random matrices. When we discuss certain large N limits in 
the GUE such as the bulk and the edge scaling limit, the nontrivial step reduces to the 
scaling limit of the Hermite polynomial in (14.20(1 . The same thing happens in our case: 
the only nontrivial step for getting (15.391) is the asymptotics of the functions fik(x\ t) (12. ip 
and 4>k(x]t) (14.ip . Based on the saddle point method, we obtain the following results of 
which the proof is given in Appendix ICl 


Lemma 18. 


'Vji N^^C(N) 

& CW) MX '' t] = (2 7t) i/ 2 Mx ‘ ; t = Aife - A), t = 1.2. 

where we set x t as (15.411) and k and t as 

AT/ 2 


(5.42) 


k = N + 


(2 7t ) 1 / 2 


A, t = v / TiV + 2 7 2 U. 


(5.43) 


The constant C(N ) is represented as C(N ) = e^ =1 Cj in terms of Ci, ■ ■ ■ , C 5 defined 
by (IC.lOp . (IG.14P and (IG.16P in Appendix [0 


On the other hand, when we take the same limit for the other representation (14.23D . 
we can apply the saddle point analysis also to the kernel (I4.24p and can get the limiting 
kernel. But since it does not correspond to the kernel (15.37P directly, we need an additional 
step to show the equivalence between the Fredholm determinant with the limiting kernel 
and that with (I5.37P (see Sec. 5.4.3 in [IU]). 

Proof of Proposition 1171 Combining the estimate (I5.42p with the simple fact 

pXi-U plT^i-s) 

_—_ i — 12 

e Xi-U _j_ e 7 T^i-s) _|_ ]_ ’ 

under (15.341) and (15.4ip . we immediately obtain the result (15.391) 


(5.44) 

□ 


6 Conclusion 

For the O’Conncll-Yor directed random polymer model, we have established the repre¬ 
sentation (12. 7p of the moment generating function for the partition function in terms of a 
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determinantal function which is regarded as a one-parameter deformation of the eigenvalue 
density function of the GUE random matrices. 

There are some special mathematical structures behind the O’Conncll-Yor model which 
play a crucial role in deriving the relation. The first one has been the determinantal 
representation (12. lip which is essentially the one with the Sklyanin measure in [57]. Next 
we have introduced another determinantal measure in enlarged degrees of freedom (12. 19[) . 
Our main theorem has been readily obtained from a simple fact about two marginals of 
this measure (Theorem [6]) . 

We can regard our approach as a generalization of the one in [78J which retains its 
determinantal structures. To see this we needed to reinterpret the dynamics in the Gelfand- 
Tsetlin cone introduced in [78] using the weight (13.33ft supported on the partially ordered 
space Vjv 113.34)1 . Our approach is a natural generalization of [78] from this viewpoint. 
It would be an interesting future problem to find a clear relation with the Macdonald 
process ra, which is another generalization of (78) . 

Applying familiar techniques in random matrix theory to the main result, we have read¬ 
ily obtained the Fredholm determinant representation of the moment generating function 
whose kernel is expressed as the biorthogonal functions both of which are simple deforma¬ 
tions of the Hcrmite polynomials. The asymptotics of the kernel under the scaling limit to 
the KPZ equation can be estimated easily by the saddle point analysis. 


A Proof of Lemma [5- 


First we give a proof of f!2.29p . For this purpose, it is sufficient to show the case of 
m — 1, x = 0, 


L dxe ^ Mx)= r ot {j}- 

Furthermore setting e^ x — y, a//3 — b, one sees that (lA.lj) is rewritten as 

POO 

/ dy h b (y ) = n cot nb, 


(A.l) 


(A.2) 


where h b (y) = y~ b ^/(y — 1) and we take the branch cut of h b (y) to be the positive real 
axis. Hence here we prove (1A.2I) . We set the contour C as depicted in Fig. [2] with a — 1. 
From the Cauchy integral theorem, we find 


>c 


dy h b (y) = 0. 


(A.3) 


Dividing the contour C into Q, i — 1, • • • , 6 as in Fig. [2j we find that by simple calculations, 


c lim _ / dyh b (y) = - lim e 2mb 

S —^0,6—^0 / s~i S —^0,e—>-0 

D ’ J ^1 D ’ 


/ dyh b (y)= / dy h b (y ), 
r c 4 Jo 
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lim / dy hb(y) = lime 


2-nib 


e—7-0 


'C 2 dy 


e —>0 


dyh b (y ) = -m, 


'C 5 


lim / dyh b (y) = lim / dy h b (y) = 0, 

Jc 3 s ^° Jc 6 

where note that the factors e 2mb s come from the cut locus of y~ b . 
From (IA.3I) . (1A.4I) . we get 


(A.4) 


fl-> oo 7=1 C '7 


e- 2mb ) / dyh b (y) - (1 + e~ 2mb )ni, 


(A.5) 


R—^oo J — 

which leads to (1A.2I) . 



Figure 2: The contour C on C, where a G (0, oo). It consists of the paths C i, • • • , Cq. 


Next we give a proof of (12.3011 . For this purpose we first show the following relation. 
Let j — 1, 2, • ■ •, x G (0, oo), be 


Then we have 



1 (log w) j 1 
x — w w + 1 




1 

X + 1 


Xj( log®), 


(A.6) 


(A.7) 


where Tk(x) (k — 0,1, 2, • • ■) is a fcth order polynomial of x where the coefficient of the 
highest degree is 1/k. This relation (IA.7[) can be derived by considering the integration 
of m,j(w,x ) := (logtc) J /((x — w)(w + 1)), x > 0 with respect to w along the contour C 
in Fig. [2] with a = x and R > 1. 
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Note that 


7 n 


IV 


/ dwmJw: x) — 2ni 
lc X +1 


(A.8) 


where RHS corresponds to the residue of rrij(w; x)atw = —1. As in the previous case (I A. 41) . 
one easily gets 


lim / dw mAw, x) — L+Ax), 

n°i^° Jc, 


lim 

<5—>0,e—>0 

Ft —>-oo 


c 4 


dwrrij(w; x) — — I dw rrij(we 2m ; x) = 


/ (log w + 2'KiV 

/ dw 7-77-— 

1 0 {x — w)(w + l) 


k =0 


- El (2 


lim [ dwmjiw] x) = m MA , lim/ mi ( TO; x)dw = iri + 


e—>0 


'C 2 


a; + 1 e->o 


C 5 


X + 1 


lim / dvj x) — lim / dwrrij(w;x) = 0. 


R—yoo 


'C 3 S ^°JCe 

Substituting (1A.9() into (1A.8I) . we hnd 

6 

'mj[w;x) 


(A.9) 


27TZ- 


7T7 




6 f 

lim y / dw: 

>0,e—>0 jf-i 


X “I - 1 —>-0,fc—fU " / s~i 

R—too k=l k 


Ij+i{x) - (fr) ( 2?ri y k/ k+i(x) + ((loga;) J + (logx + 2vri) J ) . 

k =0 b / X ' 


(A.10) 


Thus we obtain 

(logo;)- 7 + (logic + 2ni) j — 2{j\ i) j 


3-2 


tjfr) = 


2j(x + l) 


k =0 


--E 


(A.ll) 


which leads to (1A.7D . 

Here we show (12.301) . We hnd that (I2.30j) is rewritten as 

J* B {m) J F (x) = q m Jf(x), (A.12) 

where q m (x) is dehned below (12.301) and the functions J F {x) and Jb(x ) on M + are dehned 
by J F {x) = l/(ic + 1) and Jb(x) = 1 /(3x. 
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We prove (1A.12D by using (1A.7D and by mathematical induction: suppose that (IA.12|i 
holds for 77i = N — 1. Then we get 


N -2 


J *b N) Jf{x) = ± I dy ~~~ J*b N = a N( n-i\M y) + E 


If, 1 ,*(JV- 1 ) T / \ 1 r , \ °k 


P Jo x - u 


(3 N (N — 1)! 


k =0 


pk+l 


(A.13) 


where Ck(k = 0,1, ■■■, N — 2) is the coefficient of x k in q N _i(x) and in the last equality 
we used the assumption for the mathematical induction and (IA.6D . Considering (1A.7I) . we 
arrive at (IA.12p . □ 


B Proof of Lemma [TO. 


To show Lemma [Till we will use the following identity. For (aq, • • • , %_i) G R^ 1 satisfying 
X\ > • • • > x jy—i and (y i, • • • , a/jv) G R^, we have 


N 


N 


Y sg no- Y[ 1>0 (Xj-l - ya(j)) = Y sgna n 1>0 ( X 3-1 - y°U)) X > 0 (?A(J-1) - Xj- 1) • 

(TGSjv 3 =2 <t£S n j =2 

(B.l) 


where Sn is the permutation of (1,2, ••• , N). For the proof of (IB. ID . it is sufficient to 
show for m — 1, 2, • ■ ■ ,N, 


N 


Y “ W(j)) • n 1 >o(^(A:-l) - x k- 1) • 1>0 ( x m ~ Va(m)) = 0, (B.2) 

ct&Sn j =2 


k =2 


where, as in (IB. ID . we assume the condition aq > x 2 > ■ ■ ■ > aqv-i- This can easily be 
obtained by noting that 


N 


^ ^ SgllCT | 1 >0 ( x j —1 y<r(j)) ■ l>o(2/cr(A;— 1) Xk—l) • 1>0 ( x m ya(m)) 

ctGS n j =2 


k =2 


iV 


^ ^ Sgna 11 1 >o(A? — 1 y<r{j)) ’ J_ l>o(l/(T(A;—1) x k—l) ' l>o(*^m 2/<r(m)) 


o&Sn 3 =2 

j^m. 

N 


k =2 


^ ^ Sg n a 11 l>o (^--1 y&(j)) J_ l>o(?/(T(fc— 1 ) x k—l) ' l>o(*^m ?/o-(m))l>o(*^m 2/<r(m+l)) 


o-GSat i=2 

j^m, m+1 


k=2 


= o 


(B.3) 


where in the first equality we used the fact that the factor l>o(^ m -i m 'y<j(m)) can be omitted 
in this equation thanks to the factor l >0 (x m — y a (m)) with the condition Xm-\ > x m and 
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the last equality follows from the fact that the term with o cancels the term with o' where 
o' is defined in terms of o as o'[m ) = o[m + 1) and o'(m + 1) = o(m ) with o'[k) = o(k) 
for k 7^ m, m + 1. 

Using (IB.2D . we have for x\ > X 2 > ■ ■ ■ > xjv-i 


N N 

Y sgn^n 1 - (-U-I - y<r{j)) = Y sgna n 1>0 (-U-i - Vo(j)) ■ (! - !>o(^i 

o-eSjv j= 2 o-SSiv J=2 

N 

= Y sgna n 1>0 Oh ?- 1 - ?Md) • l>o(j/<r(l) - ®l) 

o-eSjv i=2 

N 

= Y s g ncr n 1>0 Oh/-i - 2/o-C#)) ' !>o(?/a(i) - a?l)(l - l>o(^2 - y ff (2))) 

cre5jv i=2 

V 2 

= x] sgncr n i>o o^-i - 2/o-c#)) ■ n i >o(^(u - ®fc)- 

ctGSn j —2 fe=l 


2/(j(l))) 


(B.4) 


where for the first and the third equality, we used (1B.2|) with m = 1 and m = 2 respectively. 
Performing the procedure in (IB. dll repeatedly, we arrive at (IB.ID . 

Now we give a proof of the lemma by the mathematical induction. The case N — 1 in 
is trivial. Suppose that it holds for N — 1. Then noticing 


N k 


N 


iyvfev) = nil 1: 


>o I 4_i 13 - x ) = 


lyv-iGUv-i) n 1; 


(V—1) JN) 


>0 X-l - X 


3 


k= 2 j =2 

we see that LHS of (13.41)1 is written as 


i =2 


(B.5) 


N 


Y Sgll(J (Ar) 1 VKr—i fev-l) II 1>0 ( X 

aU)eSj,j=l,-,N 3 =2 


(V—1) _ (N) 


N 


/ „ Sgncr J-Viv-! UZiV-lJ Sgncr XX 1 >0l a; j-l X (T (N Hj) 


V)eSj,j=i,— ,N 


3 =2 
V 


Z sgncr (iV 1} IgtGejv-i) ' sgncr (Ar) JJ 1 >0 (x^- x { ^ ){j) 


o-(j)£Sj,j=l,-,N 


3 =2 


V 


Z sgn^ (Ar) iGTfeSr-r) ]J l>o ~ X S ( J > 


rV'>eS i ,j=l,-,N 


3= 2 


(B. 6 ) 


where in the second equality we used the assumption for N — 1. Note that in the rightmost 
side of (IB.61) . the condition x^n-i)^ > x ^n-i)( 2 ) > ''' > holds f° r the support 


( 2 ) 


>(JV— 1 ) 
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of IgtGejv-i)- Thus we can apply (IB. ID to the rightmost side. We see that it becomes 


N 


ST __ _(JV) ( « \ TT 1 (AN-l) _ JN) \ 1 (N) _ (N- 1 ) 

/ , Sgncr J-GTUliV-lJ -*->0 I X a( N ){j)J i>0 \ X a( N )(j-l) X a ( N ~ 1 )(j- 1) 


aU)£Sj,j=l,~,N 


3 = 2 


sgncr (7V) 1 G t(^at)) 

o'JkS,.j I.....V 

which completes the proof of Lemma [101 


(B.7) 


□ 


C The saddle point analysis of t ) 


In this Appendix, we give a proof of (j5.42[) based on the saddle point method in a similar 
way to Sec. 5.4.3 in pLOj. Here we deal with the case of general Y while the case of Y = 0 
was considered in [TO] . We focus mainly on the limit about (12.ip in (15.42)1 since 

the case (j>k(x‘,t) (14.ip can also be estimated in a parallel way. Changing the variable as 
w = —iVNz, (I2.1jl becomes 




Vn 

27 ri 


d ze M z ^ x \ 


(C.l) 


where 


,zH 


f N (z;t,x ) = —VNzx + + {k — N) \og(\/~Nz) — N\ogT(VNz). (C.2) 


Substituting (15.41)1 and (15.431) into (1C.2I) . we arrange the first three terms in ascending 
order of powers of N as 

~ 2 t 

- VNzXi + N— + (k — N) log (VNz) 

= N 3 /' 2 log N • | + N 3 / 2 (-z + ~ \ log Tj + N 2 Yz 2 - 2 y 2 YT~ 1/2 z - 


+ N 1/2 log IV- 


X 


2(2 7 T ) 1 / 2 


+ 


n i/2 ( 7r 


A 


V u z ~ 7t6 " + 1tY ' 2z ~ ^ YZ + (2 7t )V2 


log 2 . 


Using the Stirling formula 


, rn / \ , log 27m 1 _o. 

logf (n) = nlogn — n -1-b (Jin ) 

2 12 n 


(C.3) 


(C.4) 


for the last term in (|C.2[) . we have 


— NlogT(\/Nz) = —IV 3 / 2 log N • X -b N 3 ^ 2 (,z — z log z) + ^ logN + 
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- n1/2 — + °( n ~ 1/2 ). 

Thus from (1C.311 and (1C.5(1 . /jv(z) (1C.211 can be expressed as 

f N {z- 1 , x ) = N 3 ' 2 f(z) + Ng(z) + N^ 2 h{z) + C, + 0{N~ 1 / 2 ), 


m = 

T l / 2 z 2 

2 

T x ' 2 z 

9(z) = 

2 

h(z) = 

1 

12^ 4 


z 

2 


+ 


2t|F 
T 1 / 2 ' 


7t 


+ (i-^-7ue,-r 2 )j^ (2TT)I/2 


A 


log z. 


Here Ci, which does not depend on z is written as 


(C.5) 

(C.6) 

(C.7) 

(C.8) 

(C.9) 


Ci = N 1/2 log TV . 2 ^ )1/2 + j log2vr\^V. (C.10) 

We note that /(z) above has a double saddle point 2 C = T^ 1 / 2 such that f'(z c ) = f"(z c ) = 0. 
We expand /(z), #(», h(z) around z c . Noting f"(z c ) = 27^, g'Oc) = 0, g"{z c ) = 
2 1t y - It 7 h'{z c ) = 7^/4 - 7 \Y + 7 T (A - & + >'' 2 ), we get 

N 3 ^ 2 f (z) = N 3 / 2 f{z c ) + ^f'"(z c )(z - z 3 ) 3 + 0(iV 3 / 2 (z - 2 C ) 4 ) 

= C 2 + N 3 ' 2 ^-{z - 2 C ) 3 + 0(N 3 ' 2 (z - z c ) 4 ), (C.ll) 

Ng(z) = Ng{z c ) + Ng'(z c )(z - z c ) + N^-^-(z - 2 C ) 2 + 0(iV(z - 2 C ) 3 ) 

= C 3 - N - 7 2 f) (z - z c ) 2 + 0(iV(* - 2 C ) 3 ), (C.12) 

N 1 / 2 h(z) = N 1 / 2 h(z c ) + 7V 1/2 /i'(0(z - z c ) + C 4 + 0(iV 1/2 (^ - z c ) 2 ) 

= C A + N 1 ' 2 ($■ - 7+ 7t (A - £ + H 2 )) (* - * c ) + 0(iV 1 / 2 (^ - z c ) 2 ), 

(C.13) 


where C2, C3 and C 4 are 

C 2 := N 3 ' 2 f(z c ) = N 3 / 2 T~ 1 / 2 /2, C 3 := Ng(z c ) = -TV ( 

/mi/ 2 

c 4 := iv I/2 fcfe) = -iv 1 / 2 + ( 7 Jr + 7T« - i' 2 ))r (2 ^ )1/2 

Tims from (lC.6jl - (lC.13(l and under the scaling — (z — z c ), we have 


f l + log T i y 
2 + 2^ 

, AlogT’ 1 / 2 


(C.14) 


f N {z\t,Xi) = 


It 


^ ~ _ 2 ) _ 7T ^ _ {f + ^ ~ 9 ) + Z) C i + 0 V2 ) ’ 


7r 


j=i 


(C.15) 
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(C.16) 


where C \, • • • , C 4 are defined in (1C.10|) and (1C.14|) and C 5 is 


a = —- ( f 


3V- + 

Further changing the variable — 1/2 = —iv/^/T, we obtain 

5 


t ®i) = - A)u + C7,- + O (TV 1/2 ) , 

j=i 


Hence from (1C.II) and (1C.171) . we get the limiting form of 


e ^= lCjr )T^k{xi,t) ~ — / dne^ 3+ * (ei A),; = Ai(& - A), 

27T 


(C.17) 


(C.18) 


which is nothing but (j5.42|) . 

As with (IC.ip . we rewrite 4>k{'X] t) (14.ip by the change of variable v = V~Nz, 

1 r e ~f N (z;t,x ) 

0/cO; Q = -~i f dz -7-’ (C.19) 

where /jv(z; t, x) is given in (1C.21) . Applying the same techniques as above to this equation, 
we get the result for 0 fc (a;; t). 
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